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Superradiance typically requires preparation of atoms in highly entangled multi-particle states, 
the so-called Dicke states. In this paper we discuss an alternative route where we prepare such states 
from initially uncorrelated atoms by a measurement process. By measuring higher order intensity 
intensity correlations we demonstrate that we can simulate the emission characteristics of Dicke 
superradiance by starting with atoms in the fully excited state. We describe the essence of the scheme 
by first investigating two excited atoms. Here we demonstrate how via Hanbury Brown and Twiss 
type of measurements we can produce Dicke superradiance and subradiance displayed commonly 
with two atoms in the single excited symmetric and antisymmetric Dicke states, respectively. We 
thereafter generalize the scheme to arbitrary numbers of atoms and detectors, and explain in detail 
the mechanism which leads to this result. The approach shows that Hanbury Brown and Twiss type 
intensity interference and the phenomenon of Dicke superradiance can be regarded as two sides of 
the same coin. We also present a compact result for the characteristic functional which generates 
all order intensity intensity correlations. 
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I. INTRODUCTION 

Modification of spontaneous decay is one of the fun¬ 
damental topics in quantum physics. In a pioneering 
paper, Dicke introduced in 1954 the concept of super¬ 
radiance, i.e., the coherent emission of spontaneous radi¬ 
ation [1]. The phenomenon is displayed by quantum sys¬ 
tems in particular correlated states, the so-called Dicke 
states, leading to profound modihcations of the tem¬ 
poral, directional and spectral emission characteristics 
of the ensemble compared to that of a single atom [1- 
16]. Even though a tremendous amount of literature - 
both, theoretical and experimental - has been published 
since, the physical origin of the phenomenon remained 
largely obscure. A deeper understanding has emerged 
recently in terms of quantum interferences among multi¬ 
ple path ways produced by systems in correlated states 
[17]. Based on this interpretation we investigate in the 
present paper a new aspect of superradiance, namely that 
it can be observed also with statistically independent and 
initially uncorrelated incoherent sources [18]. 

In our approach the production of atomic correlations 
and the corresponding superradiant behavior relies on 
the successive measurement of photons at particular posi¬ 
tions in the far field of the sources, such that the detection 
is unable to identify the individual photon source. In this 
case, the initially fully excited uncorrelated atomic sys¬ 
tem cascades down the ladder of symmetric Dicke states 
each time a photon is recorded. This is another exam¬ 
ple of measurement induced entanglement among parties 
which do not interact with each other and are separated 


even by macroscopic distances [19, 20]. As discussed be¬ 
low, recording m photons scattered from N > m atoms 
amounts to measuring the m-th order photon correlation 
function. Measuring this function thus allows (a) to pro¬ 
duce any desired symmetric Dicke state from statistically 
independent and initially uncorrelated incoherent atomic 
sources and (b) to observe the corresponding superradi¬ 
ant emission characteristics of the related Dicke state. 

The setup employed to display the superradiant emis¬ 
sion characteristics of initially uncorrelated incoherent 
atomic sources is similar to the one used in the celebrated 
Hanbury Brown and Twiss experiment to measure the 
angular diameter or the separation of stars [21, 22]. The 
detailed analysis shows that the two effects derive indeed 
from the same cause, namely from multi-photon inter¬ 
ferences appearing in the m-th order photon correlation 
function. In this way we show that Hanbury Brown and 
Twiss intensity interference and the phenomenon of Dicke 
superradiance can be regarded as being two sides of the 
same coin [18]. Additionally, in a sense, we provide the 
great utility of Glauber’s program [23] to extract com¬ 
plete quantum statistical information on radiation helds 
by means of higher order photon correlation measure¬ 
ments. Last but not least, we show how Dicke superra¬ 
diance can be simulated by the measurements of higher 
order intensity correlations on statistically independent 
and initially uncorrelated incoherent sources. 

The paper is organized as follows. In Sec. H we in¬ 
vestigate in detail the second order intensity correlation 
function for two independent atoms in the fully excited 
state. We show that the function displays the same emis- 
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sion characteristics as Dicke super- and subradiance for 
two atoms in the single excited symmetric and antisym¬ 
metric Dicke state, respectively. It is demonstrated that 
the detection of the first photon can be regarded as a 
projection of the initially uncorrelated system onto one 
of the maximally entangled Dicke states, and the detec¬ 
tion of the second photon as a probe of the super- and 
subradiant emission behavior of the corresponding Dicke 
state. In Sec. Ill we generalize this idea and investi¬ 
gate the TO-th order correlation function for N arbitrary 
fully excited independent atoms. Considering all possi¬ 
ble multi-photon quantum paths leading to a valid m- 
photon detection event we show that if {m— 1) detectors 
are placed at the same position we obtain again a su¬ 
perradiant intensity distribution. In Sec. IV we compare 
the results of Sec. Ill to the intensity measurement of 
N atoms being prepared in symmetric Dicke states. We 
demonstrate that the intensity distribution of a Dicke 
state with N — (m—l) excited atoms corresponds indeed 
to the m-th order correlation function for the initially 
fully excited state - what proves the state projection by 
photon subtraction. In particular for m = N, we pre¬ 
pare effectively the Dicke state with a single excitation 
and the measurement of the N-th order correlation be¬ 
comes equivalent to the measurement of for the state 
with a single excitation. Sec. V introduces the character¬ 
istic functional, which can be used for calculating higher 
order intensity correlation functions in a compact way. 
For (to — 1) detectors at the same position identical re¬ 
sults as in Sec. Ill are obtained. In Sec. VI we finally 
conclude. Note that in a subsequent paper we will apply 
the ideas of this paper to classical sources and discuss 
in detail the aspects of superradiant emission from such 
kind of emitters. 


II. HANBURY BROWN AND TWISS EFFECT 
AND DICKE SUPER- AND SUBRADIANCE 
WITH TWO RADIATING ATOMS 

The classical Hanbury Brown and Twiss effect was 
measured with classical thermal sources. The measured 
quantity was the spatial intensity intensity correlation 
(/(ri)/(r 2 )) as recorded by two detectors located at ri 
and r 2 . Now we consider the Hanbury Brown and Twiss 
measurement with two independent fully excited atoms 
at Ri and R 2 , separated by distance d much larger 
than the wavelength A of the emitted photons, such that 
the dipole-dipole coupling between the atoms can be ne¬ 
glected. With reference to Fig. 1 we consider the two 
detectors positioned in the far field in a circle around the 
sources. For simplicity we suppose that the atoms and 
the detectors are in one plane and that the atomic dipole 
moments are oriented perpendicular to this plane. The 
measured quantity is 

G(2)(ri,r2) = (£'(-)(ri).£'(-)(r2)-£’W(r2)-f'+Hri)) , 

( 1 ) 


where (ri) and (ri) are the positive and negative 
frequency parts of the electric field operator, respectively. 

Each two-level atom is represented by the spin 1 /2 op¬ 
erators (^ = 1,2). The field at the detector position 
can be related to the atomic operators via the well known 
relation [24] 




^;2 gd^ri-wt) 

ri 


= X (m X p)) 

1^1,2 


( 2 ) 

where p is the transition dipole moment, ni = ri/|ri| is 
the direction of observation, is the free field operator 
and — isy^ = Sx'^ + isy^) is the atomic 

lowering (raising) operator for atom 1. 

With the initial state |<i)) = |e, e) and using Eq. (2), the 
intensity intensity correlation function can be calculated 
to be 


G^^'>{ri,r2) = g‘^2 2 l(”i ^ P)n(n2 x p)|^G(^^(ri,r2) , 
c rir 2 

( 3 ) 

where 

G^^Hri,r 2 ) = 2[l-hcos (^(ni - n 2 ) • (Ri - R 2 ))] , (4) 

with rij = Yj/\vj\ [j = 1, 2) denoting the two unit vectors 
pointing towards the two detectors. In what follows we 
work with G^^^ which effectively is defined by 


G ( 2Hri , r2 ) = (e^-\vi)E^-\v2)E'^+\v2)E^+\v^)) , 

( 5 ) 

with 


E(-)(r)j =E(+)(r) = ^ 


(1) 


( 6 ) 



FIG. 1. Considered setup: Two identical two-level atoms, 
separated by a distance d ^ X, are placed at positions R;, 
I — 1,2; the light scattered by the sources is measured by two 
detectors, located at positions r^, j = 1,2 in the far field. 
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For brevity we will drop the bar froni Note the 

presence of fringes in the intensity intensity correlations 
(Eq. (4)) as the detector positions are varied. 

The normalized form of the intensity intensity correla¬ 
tion reads [25] 


g(2)(ri,r2) 


with 


G^^\vI,V2) 

G(i)(ri)G(i)(r2) 

^[1 -h cos (^(ni - n 2 ) • (Ri - R 2 ))] , 

( 7 ) 


G(i)(r) = <^G(-)(r)£;(+)(r)^ = 2 . (8) 


Note that there are some similarities to the Hanbury 
Brown and Twiss result for thermal sources [26] 

5^^Lmai(l’l,I-2) = l+|7(ri-r2)|2 , (9) 


where 7 (ri — r 2 ) is the complex degree of first order coher¬ 
ence between the two thermal sources. However, unlike 
Eq. (9), we could have 

—>■ 0 for single atom sources, (10) 


i.e., the probability of detecting one photon at each of 
the detectors Di and D 2 can be zero in case that 

-(ni - n2) • (Ri - R 2 ) = tt . (11) 

c 

This is the well known Hong-Ou-Mandel effect [27] with 
two identical single photons - with the single photons 
being emitted by the identically excited atoms. We 
have thus established a connection between the Hanbury 
Brown and Twiss effect and the Hong-Ou-Mandel effect 
for two identical single photon sources 

Dicke superradiance and subradiance require that the 
two atoms are prepared in an entangled state. Let us 
consider the initially prepared state 

l^) = ^(|e,5) + e*'|5,e)) , (12) 

where <5 is a parameter defining for (5 = 0 (tt) a symmetric 
(antisymmetric) state. 

From Eq. (12), the atomic expectation values are given 

by 

(4‘>4^>) = y“, (7>=“) = i. (13) 
Using Eq. (13) the intensity of the emission becomes 

= (4‘>7‘>) + (4>4>)+2»(4'4')e<f 

= 1-1- cos[d -I- —ni • (Ri - R 2 )] . 

(14) 


^ A detailed discussion of this aspect for arbitrary numbers of sin¬ 
gle photon sources will be presented in [28] 


The intensity distribution exhibits a fringe pattern as the 
detector position is scanned. In particular, for a detection 
perpendicular to the line joining the atoms, we find 1 = 2 
for (5 = 0 and / = 0 for (5 = tt (a similar result is obtained 
if the two atoms are confined to a region much smaller 
than the wavelength A, i.e., |yni • (Ri — R 2 )| « 1). It 
should be borne in mind that there is only a single ex¬ 
citation in the system. Thus for perpendicular detection 
(or for two atoms separated by d <C A), we obtain Dicke 
superradiance and subradiance for (5 = 0 and 6 = tt, re¬ 
spectively. For these values, the states of Eq. (12) are 
the Dicke states \j,m) with j = 1, m = 0 for (5 = 0 and 
j = 0, m = 0 for (5 = TT. 

If we compare Eq. (7) with Eq. (14) we see that the 
fringe patterns in Eq. (14) and Eq. (7) are identical if we 
take 


(5 = --n2-(Ri-R2) . (15) 

c 

Clearly there must then be a deep reason for this to hold. 
It implies that a measurement of G^^^(ri) with the initial 
state given by Eq. (12), along with Eq. (15), is the same 
as a measurement of G*-^^ (r 2 , ri) for the initial state |e, e). 
This clearly has an important consequence as it shows 
that the physics of the entangled state Eq. (12) can be 
explored via a measurement of G*-^^ on the state |e, e). 
Thus the measurement of G*^^^ bypasses the need for the 
preparation of the entangled state Eq. (12) as the super¬ 
radiant and subradiant characteristics of Eq (12) can be 
studied via the G^^^ measurement on the state |e, e). 

We next present a mathematical reasoning for this 
finding. Let us write G^^^(r 2 ,ri) for the initial state 
|4>) = |e, e) in the form 


Gi"^(r2,ri) = ($|F;(-)(r2)G(-)(ri)F;(+)(ri)£;(+)(r2)|$) 


(Tr[pF;(-)(ri)E;(+)(ri)])GW(r2) . 


where we introduce 

E(+)(r2)|$)($|E(-)(r2) 


P = 


G^^\r2) 


(16) 

Tr[p] = 1 . (17) 


and 


GW(r2) = ($|E(-)(r2)F;(+)(r2)|$) . (18) 

Rewriting Eq. (16) in the form 

G$^(r2,ri) = G'^^^(ri)G^^^(r2) (19) 

leads to the following important result: A measurement 
G|^^(r 2 , ri)/G^^^(r 2 ) on the state |<i)) = [e, e) is equiva¬ 
lent to a measurement of Gj;^^(ri) on the state p which 
is determined by the location of the detector D 2 . Note 
that according to Eq. (19), we can identify 

G^^^(ri) =G^^_^^(r2,ri), 


( 20 ) 
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as the conditional intensity intensity correlation, i.e., the 
joint probability of detecting one photon at ri condi¬ 
tioned on the detection of a photon at r 2 . The state p 
can be found using Eq. (6) in Eq. (17) with the result 


p= lx)(xl , 

lx) = ^ {\g, + |e, 

We have thus proved exactly the isomorphism between 
G'^^(r 2 ,ri) and G^^\ri). The explicit form of the 
normalized conditional intensity intensity correlation 
G$ c(i' 2 ,i'i) reads (cf. Eq. (20)) 


Gi^,i(r2,ri) 


1 -I- cos 



n 2 ) • (Ri 


Clearly, according to Eq. (22), we have 




which are clear signatures of superradiance and subradi¬ 
ance from the two atom system. A great advantage of 
this strategy is that we need not to prepare the system 
in an entangled state. 

The special case discussed above is applicable more 
generally, and in particular is not restricted to single 
photon sources. To demonstrate that let us think of the 
radiation field in the state p. Then Gp '{r 2 ,ri) can be 
written as 

G(,2)(r2,ri) = Tr[p£;(-)(r2)E;(-)(ri)E;(+)(ri)E;(+)(r2)] 

= G^^^(ri)GW(r2) = G^2)(r2, ri)GW(r2) , 

(24) 

with 

G(2)(r2,ri) = G;,')(ri) , (25) 

where 

, _ E(+)(r2)pA(-)(r2) 

^ Tr[pE(-)(r2)E(+)(r2)] ’ ^ ’ 

Thus the conditional intensity intensity correlation 
Gp^c(r 2 , Ti) is again the same as G^^^(ri) of the projected 
state p. Note that is the annihilation operator 

and hence p is obtained by subtracting a photon from the 
state p. Thus G^p\r 2 , ri) can be thought of providing in¬ 
formation on the photon subtracted state p. It should be 
born in mind that the projected state p depends on the 
location at which the photon is detected. 


III. m-TH ORDER CORRELATION FUNCTION 
FOR N SINGLE PHOTON EMITTERS 



FIG. 2. Considered setup: N identical two-level atoms, 
separated by a distance d ^ A, are placed along a chain 
at positions Ri, I = 1,...,A; the light scattered by the 
sources is measured by m detectors, located at positions r^, 
j = 1,..., m, in the far field. 


with upper state |e/) and ground state \gi), assuming that 
(to — 1) spontaneously emitted photons are detected at 
ri and the m-th photon at r 2 . The to detectors are sup¬ 
posed to be located at positions Vj, j = 1,... ,to, in the 
far field in a circle around the sources and the emitters 
at positions Rj, I = 1,..., A^, along a linear chain with 
equal spacing d 3> A such that the dipole-dipole cou¬ 
pling between them can be neglected (see Fig. 2). For 
simplicity we suppose again that the emitters and the 
detectors are in one plane and that the atomic dipole 
moments of the transition |e;) —>■ \gi) are oriented per¬ 
pendicular to this plane. Note that we assume a single 
photon counting regime, where all m photons have to be 
detected for a valid measurement. That is via postse¬ 
lection we keep only those events where all to detectors 
click, all other events are dropped. This is not a hand¬ 
icap as what matters for our theoretical analysis is the 
conditional measurement of the m-th photon. 

For arbitrary detector positions the TO-th order corre¬ 
lation function is defined as 

m 

:) , (27) 

i=i 


where (: ... :) denotes the (normally ordered) quantum 
mechanical expectation value. Due to the far field condi¬ 
tion, i.e., the inability to identify the particular photon 
sources, the electric held operator at rj is given by [20] 




N 

"T- 

1^1 


(28) 


In this section we generalize the ideas of Sec. II and 
derive the m-th order correlation function for N iden¬ 
tical single photon emitters (SPE), e.g., two-level atoms 


Here, S; = \gi){ei.\ is the atomic lowering operator and 


(fiij = k 


Tj ■ Ri 


= I kd sin 9j 


(29) 
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denotes the optical phase accumulated by a photon emit¬ 
ted at R/ and detected at rj relative to a photon emitted 
at the origin (cf. Fig. 2). Note that for simplicity we 
define as in the previous section the field and hence all 
correlation functions of m-th order dimensionless; the ac¬ 
tual values can be obtained by multiplying with m 
times the intensity of a single source. 

Starting with all atoms in the excited state, i.e., in 
the uncorrelated state I^at) = Hjli \^i)^ from 

Eq. (27) for the m-th order correlation function 


G, 


r 1 


N 


N 


N 


EE E H' 




\9a,) 


<7i—l a2 — l 
<71^(72 


; i<m 


N 


N 


N 


= E E ■ 

cri—1 (T 2 —<yi-\-l 


E E n< 

:(7rn-l + l j — 1 




(30) 

Here, |||t/’)|P = (V'lV') defines the norm of the state vec¬ 
tor I?/;), |...| abbreviates absolute values, and the expres¬ 
sion denotes the sum over the symmetric group 

Sm with elements cri,...,crm- In Eq- (30) the products 
represent m-photon quantum paths with 
phases accumulated by m photons emitted 

from m sources at Rcr^ and recorded by m detectors 
at Tj. Since the particular source of a recorded pho¬ 
ton is unknown we have to sum over all possible combi¬ 
nations of TO-photon quantum paths, what is expressed 
by the sum ^^ 1=1 in the second line 

of Eq. (30). Hereby, the additional condition tJi 7 ^ aj 
{i ^ j) ensures that each source emits at most one pho¬ 
ton. Considering that several combinations of m-photon 
quantum paths lead to the same final atomic state and 
thus have to be added coherently, we end up with the 
modulus square in the third line of Eq. (30). Hereby, for 
the (^) different final atomic states, the corresponding 
transition probabilities |...p have to be summed inco¬ 
herently, what results in the first combination of sums 



FIG. 3. (Color online) m-photon quantum paths of N statis¬ 
tically independent SPE where (m — 1) photons are detected 
at ri and one photon at r 2 (see Eq. (31)). The interference 
signal (ri,..., ri, r 2 ) can be reduced to a superposition of 
m indistinguishable yet different m-photon quantum paths. 
In the figure the explicit case of N = m = 4 sources and 
detectors is shown. 


of the third line of 

Eq. (30). 

We next consider that (m — 1) detectors are placed 
at the same position ri and the last detector at r 2 (see 
Fig. 3). Under these conditions we can make use of the 
fact that out of the m! m-photon quantum paths, denoted 
by n^i in Eq. (30), (m—1)! quantum 

paths are equal since (m — 1 ) detectors are placed at the 
same position. In this case Eq. (30) takes the form 


Gul (ri,-,ri,r2) 




-i 12 






— (fm — ( m{N-m) m (m-1) sin 

— [{rri 1 ).) 1 N(N-l) 


'(N'i 


{N-1) sin 2 ( '^ii-'^i 2 ) 


(31) 


with X:<ri.....a„ = l = 2Za^ = l2Za^=a^ + l 

(7i<...<<7m 


N 

■ ■ ■ 




Finally we obtain for the m-th order correlation function 
[18] 


Gi™’(ri,...,ri,r 2 ) 

_ Af!(m-1)! ( m-1 sin^(jV "^12 ) \ 

{N-m)\ 1 N-1 ' N(N-l) j I ■ 

According to Eq. (32), the m-th order correlation func¬ 
tion as a function of r 2 displays maxima and minima with 
a visibility of 


VsPE = 


m — 1 

m + 1 - ^ ’ 


(33) 


and the central maximum having an angular width of 

27r 


Nkd ’ 

where we have used the relation 


(34) 


5ipi2 = kd cos 6*2 302 = kd502 ■ (35) 

For growing numbers of emitters an increased focusing 
of the m-th photon in the direction of r 2 = ri is thus 
observed. 

Note that for m = 1 , the visibility Vspe vanishes what 
illustrates the fact that there is no preferred direction 
of emission of the first photon, i.e., the atoms emit the 
light incoherently. For m = 1 all the atoms are in the 
excited state and it is known that such a state for short 
times does not show superradiant behavior. However, for 
m > 1 , a superradiant emission and collective decay of 
the m-th photon sets in if it is preconditioned by the mea¬ 
surement of (m- 1 ) photons along a particular direction. 
For m = N a maximum visibility of Vspe = 100% is 
obtained. Note further that in the latter case the height 
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of the central maximum grows as On the other 

hand, we find that 

= )„„ = ((m - 1 )!)“ (« - m + 1 ) ^ 

(36) 

Hence the normalized m-th order correlation function 
( G^^\ has a central maximum which is pro- 

^ \ ^ / ipi2 

portional to N for m = N. The scaling of the maxi¬ 
mum as N and the angular width as 1/iV are the typ¬ 
ical features of superradiance from a Dicke state with 
one excitation. This then suggest a close connection be¬ 
tween the m-th order correlations of a fully excited state 
with the superradiant character of the suitably prepared 
states with (m — 1) deexcitations. In the next section 
we demonstrate this relationship in detail, following the 
procedure we used in Sect. II for two fully excited atoms. 
Note that temporal aspects of Dicke superradiance are 
not addressed in our analysis as they require a separate 
study where one has to solve the superradiance master 
equation numerically [29]. 


IV. SUPERRADIANCE AND m-TH ORDER 
INTENSITY CORRELATIONS 


In order to see the connection between superradiance 
and the m-th order intensity correlation function, we 
write G^™^(ri,...,ri,r 2 ) as 

Gi™^(ri,...,ri,r 2 ) 

= Tr[p(—i)G(-)(r2)G(+)(r2)]Gi™-')(ri,... .n) , 

(37) 

where we defined 






- I Ikm—l) I , 


i.e., 


(38) 


Gi7(ri,...,ri,r2) 




(39) 


This shows that the variation of G^^(ri,..., ri, r 2 ) with 
respect to r 2 , i.e., 02 , is a function of the state j'l'm-i)- 
For simplicity we can set 0i = 0. In this case we have 


F;(+)(ri) = J- , 

i 


(40) 


where J_ is the collective lowering operator as introduced 
by Dicke [1]. Using Eq. (40) in Eq. (38) we obtain 




(J)”^-i|$jv) 


1 ^ ^ t u\ 
= |-,--(m-l)) 


(41) 


where | ^, M) is the Dicke state with cooperation num¬ 
ber N/2 and M the projection of the angular momen¬ 
tum on the z-axis, i.e., J^j^, M) = (|[J+, J_])|^, M) = 
M\y,M). This means that the state of Eq. (41) rep¬ 
resents a symmetric Dicke state where (m — 1) atoms 
are in the ground state and N — (m — 1) atoms 
are in the excited state. The normalization factor 
G^~^^(ri,... ,ri) in Eq. (38) for 0^ = 0 is thus equal 
to ff,f|(J+)™-i(J_)™-i|f,f) = (J)^J((m- 1 )!) 2 , 
so that Eq. (37) can also be written as 


Gi™Hri,...,ri,r2) 

= (vI-„_i|E(-)(r 2 )E(+)(r 2 )| 4 /„_i) J ((m - 1 )!)^ 

= y ((m- 1)!)2 , 

(42) 

with 


h^.^ir2) = («'„-i|E(-)(r2)E;(+)(r2)|«'„_i) . (43) 


Note that is the intensity of the radiation produced 

by atoms prepared in the Dicke state ^ ~ (rn — 1)), 
i.e. [17], 




= {N — m + 1) 


N — m 
N -I 


m -1 sin 2 (fV^)\ 
iV(iV-l) sin2(^) ] ■ 


(44) 

We thus have the result that the superradiant charac¬ 
ter of the Dicke state ] y, y — (m — 1)) can be simulated 
by measuring the m-th order intensity correlation func¬ 
tion G|^^(ri,..., ri, r 2 ) for 0 ^ = 0 of the fully excited 
system: 


Gi™^(ri,...,ri,r2) 

fV!(m —1)! /iV — m m—1 sin^(iV^)\ 

^ (Y-m)! yN-l ^ N{N-1) J ' 

. 

The first (m—1) detections prepare the system condition¬ 
ally in the Dicke state ] y, y — (m — 1)) whose behavior is 
probed by the detection of the m-th photon. Note from 
Eq. (44) that the peak height is proportional to N for 
m = N, and the width is proportional to 1/N. These 
are the characteristic features of the Dicke superradiance 
when the Dicke state has one excitation. 

Note that by putting (m — 1) detectors at 0i = 0, we 
obtained symmetric Dicke states. However, for 0i 7 ^ 0, 
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we could get timed Dicke states as introduced by Scully 
and collaborators [6] . For example a single detection will 
produce the state 


Wn^n-i = 




(46) 


in which one atom is in the ground state with a well 
defined phase factor. 


V. COMPACT WAY FOR CALCULATING 
HIGHER ORDER INTENSITY CORRELATIONS 


In this section we finally present a compact form for 
calculating higher order intensity correlations by use of 
the characteristic functional. In this way different orders 
of the correlation functions can be obtained via simple 
differentiation. The characteristic functional C[f{-)] for 
normally ordered correlations can be defined via [30] 


C)/)-)] =<^exp|t J f*{r)E^ )(r)dV 

X exp J f{r)E^^\r)d^r 


(47) 


entiation of the characteristic functional (cf. Eq. (49)) 


G'('n)(ri^r2,... .r^) 

^ ^ df*...df^dh...dfm 


/i=...=/„=o 

/i'=---=/7.=o 


(52) 

Since all atoms are independent of each other, Eq. (49) 
can be factorized and therefore changed into a product- 
form 


N 

N 

i=i 

N N / m \ 

=n(i-i/3.f)=n > 

j=i j=i \ 1=1 j 

(53) 

where we exploited the fact that all atoms are fully ex¬ 
cited. We have thus obtained a compact form for the 
characteristic functional for the radiation from a sample 
of fully excited atoms. Again, considering the superradi¬ 
ant case by choosing 

fi = fi (i = l,2,...,m- 1) 

f - f 
Jm — J2 5 


For N sources we will use the discrete version 


E^+\ri) 


E^-Hri) 




j=i 


(48) 


we can calculate the explicit form of the correlation func¬ 
tions from Eq. (53). This procedure leads to Eq. (32). 


VI. CONCLUSION 


where cij = exp ri; • } describes the accumu¬ 

lated phase of a photon traveling from source j to detec¬ 
tor I (cf. Eq. 2). In this way we obtain 


C[f{-)] 


= { exp 




1^1 


exp 




1^1 


(49) 

where fi is the discrete version of /(r). Hereby, the sum 
over the electric field components can also be written in 
the form 


N 


N 


^//£;W(rz) , (50) 


1=1 


1=1 j=i 


i=i 


The foregoing theoretical investigations show that be¬ 
yond entangled symmetric Dicke states it is also possi¬ 
ble to employ statistically independent and initially un¬ 
correlated incoherent light sources to obtain a focused 
spatial emission pattern of the emitted radiation. For 
N initially uncorrelated two-level atoms the directional 
spontaneous emission of the m-th photon is due to pre¬ 
ceding measurements of (m — 1) photons along selected 
directions, projecting the uncorrelated atoms into Dicke 
states ~ {iTT- — 1)) with N — {m — 1) excitations 

{N > m > 1). Our work thus affirms the great impor¬ 
tance of projective measurements and demonstrates how 
higher order intensity intensity correlations could be used 
to implement different physical phenomena. 


where j3j is defined as 

m 

= ■ (51) 

In general, to derive the m-th order correlation function 
GW(ri , r 2 ,..., Tm) one can utilize the functional differ- 
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